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Abstract 

We report results on a model of two coupled oscillators that undergo periodic para- 
metric modulations with a phase difference 9. Being to a large extent analytically 
solvable, the model reveals a rich 9 dependence of the regions of parametric resonance. 
In particular, the intuitive notion that anti-phase modulations are less prone to para- 
metric resonance is confirmed for sufficiently large coupling and damping. We also 
compare our results to a recently reported mean field model of collective parametric 
instability, showing that the two-oscillator model can capture much of the qualitative 
behavior of the infinite system. 

PACS numbers: 45.05.-hx, 05.45.Xt, 05.90. -^m 



1 Introduction 

Parametric resonance is a phenomenon pervading several fields of science. It occurs when the 
modulation of a system parameter causes the system to become unstable. The literature on 
parametric resonance is enormous, so the best one can do is cite a range of different subjects 
where it may play a role. Examples include mechanical systems where such resonances 
were first identified ^, ^, ^] , elementary particles , quantum dots , astrophysics Q , 
fluid mechanics |§, plasma physics 1^], electronic networks superconducting and laser 
devices [^], biomechanics |12], and even medicine Connections with chaotic systems 



have been suggested recently |14|. The simplest and perhaps most familiar example of 
parametric resonance occurs in a harmonic oscillator whose frequency varies periodically 
with time 0, g |, |. For certain ranges of modulation parameters (frequency, amplitude) 
the oscillator is unstable while for others it is stable. Even for this seemingly simple system 
the stability boundary diagram is already quite rich and complex (see below). 

A great deal of recent work has dealt with systems of coupled oscillators - again, the 
literature in this general area is enormous. However, very little attention has focused on 
systems of coupled parametric oscillators [|^, Such coupled arrays are particularly in- 
triguing because each single oscillator alone exhibits regions of stable or unstable behavior. 
Two questions arise naturally: 1) How does coupling modify the single-oscillator stabil- 
ity boundaries? 2) Are there collective parametric instabilities of the coupled system that 
are distinct from those experienced by single oscillators? In a recent paper Bena and Van 
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den Broeck [15] addressed these questions for an infinite set of globally mean- field coupled 
harmonic oscillators subjected to time-periodic block pulses with quenched uniformly dis- 
tributed random phases ("quenched" in this context means that the phase of frequency 
modulation of each oscillator is set at time t = and then remains unchanged). Within 
their mean field treatment they are able to deal with both of the questions posed above. In 
particular, they find wide ranges of parameter values that lead to collective instabilities. 

Our interest in this paper is to identify generic coupling effects that are not specifically 
a consequence of the mean field analysis, and in this context to understand how a coupled 
system of (very) few oscillators with short-ranged interactions might carry in it the seeds of 
the infinite/infinitely cross-coupled array. In particular, we seek the seeds of the collective 
parametric instabilities. We do this by studying a model of two coupled oscillators. The 
control parameter in our study is the phase difference in the periodic modulation of the 
frequencies of the two oscillators. 

In Section ^ we introduce the model and its partially analytic solution. Section ^ is a 
presentation of results as a large collection of stability boundary diagrams that convey the 
way in which stability boundaries shift, appear, and disappear with parameter changes. In 
Section Q we analyze the ways in which the seeds of the collective parametric instabilities 
found in the mean field model already appear in the two-oscillator system. A brief summary 
of our conclusions is reiterated in Section ^. 

2 The model 

Our system consists of two linearly coupled parametric oscillators whose equations of motion 
are 

xi = -uj'^[l + (l)i{t)]xi - k{xi - X2) - jxi 

X2 = -UJ^ll + (I)2it)]x2 - k{x2 - Xi) - -fX2 . (1) 

Here ojq is the natural frequency of each (uncoupled) oscillator, k is the coupling constant 
between them and 7 is the damping coefficient. The parametric modulations 0i(t) and 
02(i) are square waves with period T = 27r/c<jp, amplitude A and phase difference 9: 

(j)i{t) = Asgn{sm{ujpt)) 

Mt) = Asgn{sm{ujpt + e)) . (2) 

With A = one is left with coupled ordinary damped harmonic oscillators whose total 
energy decays exponentially to zero. In the presence of the parametric modulation, however, 
energy is periodically pumped into the system, which may or may not lead to parametric 
resonance, i.e., to an infinite growth of the amplitudes of the oscillators. Our goal is to 
determine the boundaries between these two behaviors, which will be referred to as unstable 
and stable. As a point of reference we point the reader to the appendix where we review 
the results for a single parametric oscillator. 

It is worth noting that a simple rescaling of time to dimensionless units, t' = tujQ, shows 
that Eqs. (|l]) are governed by the dimensionless parameter combinations r = uJo/iOp, A, 
k/ujQ, 7/wo and 6. Moreover, the behavior of the system is invariant with respect to the 
replacement of 9 hy 9' = 2tt — 6 (reflection around vr) since this just amounts to an exchange 
of indices between the oscillators. 
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2.1 Floquet theory 



The linearity of the equations ahows one to make use of Floquet theory to solve the prob- 
lem [^. Defining X = (xi, xi, X2, where the superscript T denotes the transpose (i.e., 
X is a column matrix), one can rewrite Eqs. (|^) in matrix form as X{t) = D{t)X{t) with a 
matrix D satisfying D{t) = D{t+T). Given a solution X{t) to Eqs. (||), the time periodicity 
of the parametric modulations thus implies that X{t + T) is also a solution. Therefore one 
can find a matrix F such that X{t + T) = F{T)X(t) and hence more generally through a 
repetition of this solution X{t + riT) = F'^{T)X{t). The long time behavior of the system 
is thus clearly determined by the eigenvalues {Aj} of the Floquet operator F, which propa- 
gates the system in phase space for one period of the modulation. The eigenvectors Vj of 
F satisfy \vj{t + nT)\ = \Xj\'^\vj{t)\, so in the limit n ^ 00 parametric resonance occurs if 
maxj{|Aj|} > 1. 

The model is appealing because for any number of oscillators (a single oscillator, or 
the coupled oscillators considered here, or even the mean field version of an infinitely 
cross-coupled infinite chain |l5|) the piecewise constant parametric modulation leads to 
a piecewise linear system whose piecewise solution is known analytically. One can therefore 
construct the Floquet operator explicitly by simply multiplying together piecewise linear 
Floquet operators. This provides an immense reduction in computational effort by skipping 
what is usually the most time-consuming task in obtaining the regions of parametric res- 
onance, namely, the numerical evaluation of the Floquet operator itself. For an isolated 
oscillator this procedure is well known |2|, |^ and reviewed in the appendix. 

For the sake of clarity, let us first analyze the frictionless case 7 = 0. Each linear interval 
is characterized by one of the four possible states of the modulations (^1, ^2) = +A), 
{+A, —A), {—A, +A), or {—A, —A). For each of these states, the solutions are of the form 
Xj{t) = AjC^^^ + Bje^^*, where the eigenfrequencies P and M follow directly from the 
diagonalization of D: 

2 _ +u;i + v/K-c^i)^+4fc2 
2 



(wf _ ^2)2 + 4^2 



(3) 



where 2 = "^ofl + 4>i,2{'t)] + k. Denoting the current state of the modulation by indices 

+ and — , notice that while Mj^j^ and = P \. are always real, M becomes 

imaginary when A> 1 while P and M_| = M |_ become imaginary when A > l + 2A;/u;g. 

Whether the system periodically alternates between harmonic behavior and one or more 
saddle nodes will thus depend on the parameter region and the phase difference 6. 

During a time interval r with fixed (/>!, 02 we can relate X{t + T) to X[t) as X(t + r) = 
f(T)X(t), where / is the piecewise linear Floquet operator: 



/ -miCp+piCm 

miP'^Sp — piM^s 

k[—Cp + Cm] 
\ k[P^Sp-M^Sm] 



-miSp +PlSjn 
-miCp+piCm 
^[~Sp + Sm] 
k[—Cp + Cm] 



k[ Cp -\- Cm] k\ Sp -\- Sjri\ \ 



k[PHp - M^Sm] 
-m2Cp + P2Cm 

m2P^ 



•Sp-p2M'^Sm 



k\ Cp + Cm\ 
-m2Sp + P2Sm 
-m2Cp+P2Cm 

(4) 
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where we make use of the shorthand notation Cp = cos(Pt), Sp = P~^sin(Pt), Cm = 
cos(Mt), Sm = sm{Mt), mi, 2 = - 2 ^'"-'^ Pi,2 = P"^ -^1,2- The Floquet operator 
is then finahy obtained as the product of piecewise Floquet operators with arguments that 
depend on the phase difference between the modulations: 



1 

vr 



1 

vr 



(5) 



The stability properties of the coupled system are determined by the magnitudes of the 
four eigenvalues Ai,...,A4 of the Floquet operator. Note that if = or = vr, the 
above expression reduces to the product of only two matrices (as it should) since /(O) 
is just the identity matrix. The addition of damping can be dealt with by noticing that 
if X{t) is a solution for 7 = 0, then Y{t) = e~'^*/^X(t) is a solution for 7 7^ 0, where 
X{t) = ^WLa 2^(^2^-72/4. 



3 Results 

In the absence of damping, Liouville's theorem insures that det / = 1, which is confirmed 
by cumbersome calculations and in turn implies A1A2A3A4 = deti^ = 1. In contrast with 
the single oscillator (see appendix), here this is not a sufficient condition to guarantee that 
the parametric resonance bifurcations occur at Aj = ±1 ||2|, Q: for coupled oscillators the 
largest eigenvalue can cross the unit circle |Aj| = 1 and hence enter a region of instability 
(parametric resonance) in other directions in the complex plane. The fourth order charac- 
teristic polynomial needs to be solved numerically to determine the instability boundaries; 
this is a computationally inexpensive procedure since we have an analytic expression for 
the Floquet operator. 



3.1 In-phase modulations 

The case = corresponding to in-phase modulation </>2(t) = (t>i{t) simplifies sufficiently 
to yield analytic results for the stability boundaries. This is the simplest situation to be 
studied, since it is possible to reduce the original four-dimensional problem in this case 
to two known two-dimensional systems. This decoupling is accomplished by changing the 
coordinate system to the reference frame of the center of mass. Defining x = (xi -|- X2)/2 
and p = (xi — X2)/2, one obtains from Eq. (|^) 

X = -Wo[l + '/'i(t)]x - 7X (6) 

p = -{ul[l + 4>i{t)]+2k]p-^p. (7) 

Each of these equations is precisely that of a single parametric oscillator, one with frequency 
ijOq and modulation amplitude A, the other with rescaled parameters 

uj'q = ujI + 2k 



A' = A(woM)'- (8) 

We can thus use the known results for the single oscillator [Q, , for which a closed expression 
exists for the boundary of the instability regions in the (r, ^4) plane (see the appendix). 
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Figure 1: Instability boundaries for 6 = in the {r,A) plane: 7 = in the first column 
while 7/0^0 = 0.3 in the second column. k/ujQ = 0.02 in the first row, while k/ujQ = 3 in 
the second row. The dashed line in the first panel is the result for A; = (single uncoupled 
oscillator). The dark points in this panel are parameter pairs to be considered in more 
detail subsequently. 



For an undamped parametric oscillator with parameters (r, A) parametric instabilities 
occur inside "tongues" that open up from integer and half-integer values of r as discussed 
in the appendix and shown by the dashed curves in the first panel of Fig. |l|. For the 
coupled problem in which we have transformed the problem to two independent parametric 
oscillators with different parameters, the instability regions are given by the overlap of the 
sets of tongues arising from each independent oscillator, one set emerging from integer and 
half- integer values of r ( "r-instability regions" ) , and the other from 

r'{k) = rJl + 2k/u:l = (9) 

Up 

( "r'-instability regions" ) . When the coupling k between the oscillators is small then the two 
sets of tongues almost overlap and one obtains in the undamped case the solid curves in the 
first panel, which show the instability boundaries when k/ujQ = 0.02. Still in the undamped 
problem, when the coupling k is large the two sets of tongues occur at different parameter 
scales. The boundary diagram for k/uj^ = 3 is shown in lower left panel of the figure. In 
either case, since now we have two independent sets of instability regions the effect of the 
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coupling k for 6 = has been to enlarge the parametric resonance regime relative to that 
of two uncoupled oscillators. 

This last result is completely general: the effect of coupling for in-phase parametric 
modulations is to enlarge the parametric resonance regime relative to the uncoupled case 
regardless of the form of the periodic modulation. The details of the phase boundaries will 
of course depend on the specific modulation function. 

Damping, even in an isolated parametric oscillator, destroys much of the intricate bound- 
ary structure and increases the regions of stability. This is also the case for coupled os- 
cillators, where increasing the damping tends to blur out the effects of coupling. These 
behaviors can be seen in the right column of Fig. |l]. The regions of instability are now 
restricted to larger amplitudes A. Note the similarity between these two figures, which 
involve different couplings but now with substantial damping, ^/ujq = 0.3. Note also that 
in the lower panel of this column (large coupling) the very first instability wedge is an r'- 
instability region while the other portions of the diagram (as well as the unstable regions 
in the upper weak-coupling case) include both r and r' instabilities. 

For the particular case 9 = these results allow us to say something about the interesting 
problem of asymptotic synchronization. In regions where the center of mass coordinate x 
is unstable but the relative coordinate p is stable (r-instability regions that do not overlap 
with r'-instability regions) the coupled oscillators are synchronized if 7 ^ (xi = X2), 
that is, the two oscillators move together about the origin with ever increasing amplitude. 
Conversely, if x is stable but p is unstable (r'-instability regions that do not overlap with r- 
instability regions), with 7 7^ the oscillators become "antisynchronized" {xi = —X2), that 
is, the two oscillators oscillate with ever increasing amplitude but in opposite directions, 
crossing one another each time they pass through the origin. Antisynchronization becomes 
more difficult to achieve with increasing coupling. If 7 = the strict equalities xi = X2 
or xi = —X2 no longer hold in the non-overlapping regions, but the difference between xi 
and X2 or —X2 is oscillatory and remains bounded. A simultaneous instability of both x 
and p involves an unstable center of mass coordinate and unbounded oscillations of each 
oscillator about this unbounded mean, which in turn involves a more complicated phase 
relation between the motions of the two oscillators. We return to this issue later. 

3.2 Out-of-phase modulations 

A wealth of very intricate results arises when 9^0. In contrast with the 9 = case, it is 
now no longer clear how to break down the problem into simpler independent components 
(even for 9 = ir) and, in particular, there is no longer a transparent way to relate the results 
to those of single parametric oscillators. 

Perhaps the most striking feature of the ^-dependence is its sensitivity: the unstable 
regions in the (r. A) plane can change very abruptly (yet continuously) with the phase 
difference. As one might expect, this sensitivity is primarily modulated by the coupling k, 
as exemplified in Fig. |2| in the absence of damping. The left and right columns of Fig. ^ 
show the instability regions in the (r. A) plane for k/ujQ = 0.05 and k/ujQ = 0.5, respectively. 
The evolution of the boundaries is shown as 9 increases vertically along the columns (in 
order to save space, this and subsequent figures consistently omit labels, using the same 
scale as Fig. || in the (r, A) plane, with the origin at the bottom left corner and tic marks 
separated by half a unit on both axes). While the most visible effects of the phase difference 
seem to concentrate on the low A region for the smaller coupling (left column), the larger 
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Figure 2: (r, A) plane as in Fig. || (origin is at bottom left corner; tic marks are separated 
by half a unit on each axis). Left column: k/ujQ = 0.05. Right column: k/u)Q = 0.5. 7 = 
for all panels. 

coupling (right column) induces a richer behavior, with (in)stabilities arising also for larger 
values oi A as 6 changes. The phase difference can therefore create regions of stability and 
instability which are absent for in-phase modulations. 

In order to address the effect of damping in this scenario, we have chosen to fix the 
coupling k and tile several (r, A) instability plots in such a way that 6 varies in the vertical 
direction while 7 varies in the horizontal direction. Fig. ^ shows such a panel for a relatively 
small coupling k/tjQ = 0.12. One notices that the 6 dependence is smoothed out as damping 



7 



Y/COo = Y/C0o = 0.01 Y/C0o=0.07 y/coo = 0.3 




Figure 3: {r,A) plane as in Fig. || for k/ojQ = 0.12 and several combinations of 9 and 7 
(origin is at bottom left corner; tic marks are separated by half a unit on each axis). 

increases, so that almost no structure is visible along the rightmost column {'j/uJo = 0.3). 
Fig. ^ shows the results for a larger coupling k/ojQ = 0.6. Note that the same tendency 
is observed: the instability boundaries become gradually less sensitive to as 7 increases. 
However, the rightmost column of Fig. ^ now shows much more structure than that of Fig. |^, 
indicating that for the same value of 7, the (r. A) plane with larger k shows a stronger 9 
dependence. 

Figs. ^ and Q suggest that 9 = tt gives rise to particularly stable behavior. This is 
in agreement with the vague intuitive notion that anti-phase modulations should be less 
prone to resonance than in-phase ones. In order to verify the extent to which this intuition 
is correct, a complementary view of these phenomena can be obtained by projecting the 
instability regions in the {k,9) plane instead. This allows us to start with uncoupled oscil- 
lators {k = 0) and observe how a given (in)stability evolves as k and 9 change. Indeed, it 
turns out to be possible to understand much of the behavior of the coupled system in terms 
of the behavior of the uncoupled system. To produce the representative results shown in 
Figs. |5| and we fix several (r. A) points in the single oscillator stability boundary diagram 
shown as black circles in the first panel of Fig. Q and study the way in which variations in 
k and 9 affect these particular states. Two of the points in Fig. ||, {r,A) = (0.8, 1.8) and 
(1.25,0.6) are stable states for the single oscillator (the first black dot touches the stability 
boundary in the figure only because it has been drawn large enough to render it visible; 
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y/COo = y/(0o = 0.01 7/0)0 = 0.07 7/033= 0.3 




Figure 4: (r, A) plane as in Fig. |T|for Aj/wq = 0.6 and several combinations of 9 and 7 (origin 
is at bottom left corner; tic marks are separated by half a unit on each axis). 

the point is well within the stable region). The other two points, {r,A) = (1.9.0.9) and 
(2.25, 1.3) lead to unstable behavior of the single oscillator. 

The first thing to be noted in Fig. || is that the horizontal axis has been rescaled in order 
to reveal the relevance of the variable r' (see Eq. (P)). The top panels focus on r = 1.9 
and A = 0.9, which is in the parametrically resonant regime for the uncoupled system. 
The small-r' portion of the figure thus represents an unstable regime. One notes that for 
r'{k) slightly above 2 the system becomes stable within a 6 interval centered around vr. 
Increasing k a little further, this stability region then evolves in a very complex pattern, 
which includes reentrant "holes" of instability. After a distinguishable gap of instability, a 
somewhat simpler band of stability arises around r' ~ 3.5 starting at ^ = vr. This band 
continues to larger values of k, with its outermost boundaries presenting a relatively simpler 
envelope than the low-/c pattern. The interesting point to be emphasized is that the band 
is perforated by gaps of instability most of which are centered precisely at integer and half- 
integer values of r'. This result is perhaps anticipated by the fact that the frequency cuq 
appears as the effective average diagonal frequency in the mean field equation of motion 
in |l^]). The gaps are eventually closed by increasing the damping (top right panel), which 
also broadens the stability band and simplifies its dependence on the phase difference. The 
bottom panels in Fig. ^ are for r = 2.25 and A = 1.3, which again is in the resonant regime 
of the single oscillator. One notices the same pattern: for lower values of k, a complex shape 
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emerges in the 9 dependence of the stable regions. For sufficiently large values of k, a band 
of stability arises which has instability gaps basically centered at integer and half-integer 
values of r' . Damping (bottom right panel) causes gaps to disappear, creating a uniform 
region of stability centered around 6 = n. 





Figure 5: Stability regions in the {k, 9) plane in the absence (left column, 7 = 0) or presence 
(right column, 'y/ujQ = 0.05) of damping. Top panels: r = 1.9 and A = 0.9. Bottom panels: 
r = 2.25 and A = 1.3. 



Fig. ^ presents what may be regarded as the opposite situation, namely, when the original 
uncoupled system is stable. The top panels show the results for r = 1.25 and A = 0.6. Notice 
that the behavior is much simpler in this case, with the original stability being disturbed 
mostly around integer and half- integer values of r' in the absence of damping (top left), 
with a relatively weak 9 dependence. The effect of damping (top right) is to suppress most 
of these instability regions, yielding a predominantly stable {k, 9) plane. The bottom panels 
show results for an interesting intermediate situation: even though the uncoupled system is 
stable for r = 0.8 and A = 1.8, this point lies in a narrow corridor between two instability 
regions in the {r,A) plane (see Fig. Q). One would therefore expect instabilities to arise 
more easily, and the immediate question is how the resulting diagram might reconcile the 
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Figure 6: Instability regions in the (fc, 9) plane in the absence (left column, 7 = 0) or 
presence (right column, 7/wo = 0.05) of damping. Top panels: r = 1.25 and A = 0.6. 
Bottom panels: r = 0.8 and A = 1.8. 



structures observed in Fig. ^ and the top panels of Fig. ^ The answer lies in a very rich 
structure in the {k,0) plane (bottom left panel): initially, very small coupling induces an 
instability for all 9. The now predominantly unstable system evolves in a manner similar to 
those of Fig. ^, a stability region being created around 9 = tt as k increases, with string-like 
gaps of instability around the usual values of r' . The difference is that there is now also a 
second band of stability centered around 9 = 0, presenting gaps at the same r' positions. 
For k sufficiently large, those two stability bands merge and one is left with a structure 
similar to that of the top panels: a predominantly stable region permeated by strings of 
instability. The bottom right panel of Fig. ^ shows the effect of damping, which significantly 
reduces the instability gaps thus greatly simplifying the picture. 

Therefore the notion that modulations operating with a phase difference 9 = ir are less 
prone to parametric resonance is thus seen to be essentially correct. Framing this statement 
more carefully, our results show that for given r and A, sufficiently large values of /j/wq and 
'y/ujQ are able to induce, in the {k, 9) plane, a band of stability centered around 9 



TT even 
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if the uncoupled oscillators are individually unstable. The width of this band can eventually 
comprise the whole 27r interval if the uncoupled system is originally stable. 



4 Collective parametric instability 



Bena and Van den Broeck ||15[ studied the stability boundaries of N parametrically modu- 
lated oscillators {xi} each coupled to all the others by the same coupling constant 2k/N (in 
our notation). The phases {9i} are initially chosen at random from a uniform distribution 
in the interval [0, 27r], remaining quenched thereafter. With square-block modulation the 
system is exactly solvable in the limit N ^ oo, where the mean- field solution becomes 
exact. The mean field equation is 

X = -Lo^[l + (l)e{t)]x -jx- 2k{x - (x)) (10) 

where x is the displacement of any oscillator in the chain, (j)o{t) is the periodic modulation 
with phase 0, and (x) = J2i^iXi is the mean displacement to be determined self- 
consistently. Bena and Van den Broeck note that the exact solution of this equation is 




+ 2.G,(t)./^'drG,(r)-.( . (11 



The propagator G{t) is known explicitly. Indeed, at t = T it is G(T) = e "''^ F{T) where F 
is precisely the single-oscillator Floquet operator given in Eqs. dAj) and (|A.10|) , but now 
with the frequencies shifted by the coupling constant 



= V^o(l±^) + 2^-774. (12) 

Note that this is exactly the same as the propagator associated with the p variable of Eq. ^ 
in the two-oscillator in-phase modulation problem, that is, the propagator associated with 
a single oscillator of frequency Wq- 

Bena and Van den Broeck identify two sorts of instabilities. One, which they call the 
"usual parametric resonance," arises from the divergence associated with eigenvalues of 
G of magnitude greater than unity, that is, with the unbounded growth of the first term 
m Eq. (0) which in turn signals the unbounded growth of the amplitude of any typical 
oscillator in the chain. The stability boundaries associated with this type of instability 
are given precisely by Eq. ( A.12| ) and are shown for the parameter choices indicated in the 



caption as the dotted curves in the top row panels of Fig. ^. In our in-phase two-oscillator 
parlance these are exactly the boundaries of the "r'-instability regions" defined in terms 
of the shifted frequency ujq (cf. Eq. (^)). The "usual" regions shrink in width and move 
toward lower r and larger A with increasing coupling fc, a behavior already exhibited in the 
context of the in-phase two-oscillator results of Fig. ||. Indeed, this instability is beyond the 
scale of the figures in the large-coupling bottom row panels. The other type of instability, 
which they call a "collective instability," is associated with the divergence of the mean (x) 
and hence of the second term in Eq. (O) . The collective instability boundaries are shown 
as dashed curves in all panels of Fig. 0. Note that the two types of instabilities may occur 
simultaneously, as seen in the instability region overlap evident in the top row panels of the 
figure. We return to this point below. Note also that with increasing coupling the system 
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Figure 7: Instability regions in the {r,A) plane for j/uq = 0.4. According to the mean 
field model [15|, usual parametric resonance occurs inside the dotted boundaries, while the 
boundaries of collective instability are depicted by dashed lines, fc/wg — ^ for the top 
panels, and k/ujQ = 20 for the bottom panels. The thin solid lines correspond to = 
(left column); the thick solid lines (right column) correspond to = vr in the two-oscillator 
model. The dark points in the upper right panel are parameter pairs considered in more 
detail subsequently. 



becomes increasingly stable, as one might expect, and that the unstable behavior becomes 
primarily collective. 

We wish to explore whether our two-oscillator model is able to capture at least some 
of the behavior of the infinite system. In particular, we would like to investigate whether 
features of the collective instabilities of the mean field model are apparent in a system of 
only two oscillators with a fixed relative modulation phase. To make the comparison we 
also show in Fig. ^ the stability boundary for a single uncoupled oscillator (thin solid lines 
in left column) and for coupled oscillators with relative modulation phase 9 = tt. 

We make the following assertions: the two-oscillator system with any value of 9 captures 
features of the overlap region of "usual" and "collective" instabilities. The purely "usual" 
regions are captured most accurately by the 9 = system, and the purely "collective" 
regime is best captured by the two-oscillator system with 9 = tt. It is therefore this latter 
system that most fully captures (with unexpected detail) the principal features of collective 
behavior of the mean field model, and does so with increasing accuracy as the coupling 
between oscillators increases. We support these assertions, particularly the last one which 
is the one of most interest to us, with the results shown in Figs. and |^. 

Clearly, the 9 = system captures the full "usual" instability regime exactly since, 
as already stated, the "usual" instability is exactly the same as the "r'-instability." This 
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identity is not restricted to the square wave modulation but holds for any periodic modu- 
lation. However, the 9 = system does not capture the "collective" instability since the 
"r-instability" condition is that associated with a single parametric oscillator with the un- 
shifted frequency ujq. Thus, for example, in the first panel of Fig. ^ the ^ = instability 
boundaries can be constructed from the combination of the dotted regimes and the thin 
solid lines (compare with the right lower panel of Fig. ||), whereas those of the mean field 
system include the same dotted regimes but now the very different dashed regions. The 
left lower panel shows an even greater difference between the 9 = two-oscillator model 
(whose "r-instability" boundaries are independent of k) and the mean field model (where 
the boundaries of the collective instability are sensitively dependent on k). Two coupled 
parametric oscillators with relative modulation phase 9 = therefore do not capture the 
collective features of the mean field model. 

Our assertion that the two-oscillator system with any 9 contains elements of the "usual" 
instabilities in the mean field model is simply a restatement of our earlier observation that 
"r'-instabilities" continue to appear even when one moves away from ^ = and all the way 
to 6* = vr. 

Consider now the two coupled oscillators with 9 = n. The stability boundaries are 
shown by the thick solid lines in the panels in the right column of Fig. ^. In the upper panel 
we observe that the first tongue approximates the region of "usual" mean field instability 
(in the non-overlapping regime) and that the remainder captures the collective instability 
boundary features surprisingly well (although it does miss the gap). We particularly point 
to the excellent fit of the leftmost boundary of this region. The agreement between the two 
models is even more dramatic in the lower panel, which corresponds to stronger coupling k. 
Again, the details are surprisingly well matched and the leftmost boundary of the region is 
captured essentially exactly. 

To further support our analysis, and to gain a clearer understanding of the difference 
between "usual" and "collective" instabilities (which are both evidently already present in 
our two-oscillator system although the notion of a collective effect is not obvious in such 
a small system), we consider the motions that might characterize the instabilities. In the 
mean field system we conjecture that in the non-overlapping "usual" instability regime 
the mean is zero, {x) = 0, but each oscillator oscillates about zero with ever increasing 
amplitude. The motion in the non-overlapping "collective" instability regime may involve 
an ever increasing mean with each oscillator oscillating about this moving mean with finite 
amplitude. This description is in accord with that of Bena and Van den Broeck [15|. The 
overlap regions may involve both an increasing mean and oscillations of ever increasing 
amplitude about the moving mean. We have not ascertained these conjectures in the mean 
field system, but present results for the two-oscillator system that support this description. 

Figure ^ shows trajectories for the two-oscillator 9 = tt system at the two points marked 
on the upper right panel of Fig. |7[ The trajectories shown are those of each of the two 
oscillators as well as the mean trajectory. The upper panel is for parameter values in the 
unstable region that is not in the "collective" regime. It is tempting to associate this with 
the non-overlapping "usual" instability of the mean field model, an association that requires 
some caution. The figure indicates that not only does each oscillator and also the mean 
oscillate about zero, but all the trajectories, including the mean, appear to diverge. This 
behavior is that envisioned in our earlier discussion of the 9 = two-oscillator system 
in the regime where "r-instabilities" and "r'-instabilities" overlap, and is an indication 
that features of both kinds of instabilities persist even at 9 = tt. We conjecture that "r- 
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Figure 8: Trajectories for the two-oscillator 6 = tt system for the parameters corresponding 
to the dark points in the upper right panel of Fig. 0: j/loq = 0.4, A = 3, and k/iOQ = A. r = 
0.25 in the upper panel and 0.4 in the lower panel. The dashed curves are the displacement 
xi of one oscillator vs time, the dotted curves are the displacement X2 of the other oscillator, 
and the thick solid curves represent the mean displacement x = (xi + X2)/2. 

instabilities" are finite size divergences not present in the mean field model (and hence we 
would not expect to see this particular type of trajectory there). In the mean field limit 
only the "r'-instability" contributions in this regime persist, becoming the non-overlapping 
"usual" instability portion of the phase boundary portrait. 

The lower panel of Fig. ^ is for parameter values in the non-overlapping "collective" 
regime. The mean indeed increases, and each oscillator oscillates about this increasing 
mean. It is interesting that the main feature of the "collective" instability, namely, an 
increasing mean displacement with individual oscillators oscillating about this mean, is 
already clearly captured by the two-oscillator model. 

5 Conclusions 

We have presented results for a model of two coupled parametric oscillators that can be 
solved for any value of the phase difference between the modulations. The other parame- 
ters of the problem are the oscillator frequency ujq, the modulation amplitude A and period 
T = 2Tr/L0p, the coupling k and the damping 7. 

We found that in-phase {6 = 0) coupled oscillator motion can be separated into two 
independent contributions. One involves the two oscillators moving together about zero ("r- 
instabilities" ) . The instability boundaries for this motion are identical to those of a single 
parametric oscillator of frequency cjq ^-nd are independent of coupling k since the spring 



15 



connecting the oscillators is never disturbed. The other involves the two oscillators moving 
about zero but exactly in anti-phase with one another ( "r'-instabilities" ) . The stability 
boundaries are sensitive to k for these motions: the system becomes more stable with 
increasing coupling. We noted that these latter instabilities are exactly those identified as 
"usual" instabilities in the mean field model [15| and that they contribute to the instability 
boundaries in our two-oscillator model for any ^, not just for ^ = 0. We also showed that 
damping shrinks the instability regimes and smoothes the stability boundaries. 

We showed that a change in 9 can substantially modify the regions of parametric insta- 
bility as projected in the usual (r = ujo/ujp,A) plane, and that these changes are strongly 
affected by the coupling between the oscillators. In general, increasing up to vr provides 
greater stability but also leads to more intricate stability boundaries. An increase in k 
and/or 7 also leads to increased stability. Alternatively, it is possible to understand the 
coupled system in terms of the uncoupled system by projecting the instability regions onto 
the (A:, 6) plane. This rendition is particularly useful to show that vr-centered bands of sta- 
bility arise for sufficiently large k and 7. We have thus identified all the trends of behavior 
in the two-oscillator model as each of the parameters is varied. 

Our most interesting insights arise from a comparison of the two-oscillator results with 
the mean field model [15| consisting of — > 00 parametrically modulated oscillators all 
coupled to one another. In this latter system two types of instabilities have been identified: 
"usual" and "collective." The "usual" instabilities are exactly our "r'-instabilities," that is, 
the stability boundaries are identical. The interesting result is that the "collective" insta- 
bilities already appear in the two-oscillator model with 6 = ir. This statement is based on 
the great similarity of the stability boundaries of the mean field and two-oscillator systems, 
especially with increasing coupling, and the specific features of the oscillator trajectories 
that typify the motions in each of these unstable regimes. It is perhaps surprising that a 
two-oscillator model can capture so much of the mean field collective behavior, and suggests 
that collective resonance in the latter may be dominated by phases quenched around tt. 



APPENDIX A - Parametric Resonance for a Single Oscillator 

A brief review of parametric resonance for a single oscillator provides results called upon in 
the body of the paper for the coupled system [§, The equation of motion is 

x + a;^(t)x + 7x = . (A.l) 

The frequency uj{t) is a periodic function with period T = 27r/wp. The simplest model of 
such a parametric modulation is perhaps the square wave 

a;2(t) =ujI[1 + Asgn{sm{ujpt))] . (A.2) 

In the undamped case (7 = 0) the oscillator frequency thus switches periodically between 
the two values ui± = wqVI i A. If ^ < 1, the system jumps between two oscillatory regimes, 
each of which is stable. If yl > 1, then uj- becomes purely imaginary, corresponding to a 
saddle node behaviour in the phase space (x,x). In this case, the system jumps between 
an unstable saddle and a stable oscillatory regime. As we shall see below, both instability 
and stability can occur for either A < 1 or ^ > 1. With damping, the relevant frequencies 
are shifted (see below). The stability of the system depends on the interplay between the 
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amplitude A of the modulation and its frequency ujp as well as the system parameters ujq 
and 7. 

The linearity of the equations leads naturally to the application of Floquet theory for 
the solution of the problem [H] . A change of variables to 



yields the undamped equation 
where 



4 

For the square wave modulation we then have 



y{t) = e^'/\{t) (A.3) 



y + Cb\t)y = (A.4) 



Cj'^it) =Lo'^{t) (A.5) 



oj = { (A.6) 

Due to the periodicity of uj{t) the two independent solutions of Eq. ( |A.4| ) must also have 
period T. One can write for the original variable 

X(t + T) = e-^^F{T)X{t) , (A.7) 

where 

and F is the so-called Floquet operator that (aside from damping) propagates the system 
in phase space for one period: 

F(r) = /_(r/2)/+(r/2) . (a.9) 

where 

(cosL^-l-i + 2^ sincj±t ^^sincZiit \ 

^s'moj±t cos oj±t — 2^ sin uj±t J 

Equation ( A.7| ) implies that 

X(t + nT) = e-'^T^F"(r)X(t) , (A.ll) 

so the long term behavior of the system is governed by the eigenvalues of e~"''^ F. If any 
eigenvalue is larger than 1 in absolute value, there is parametric resonance and the system 
is unstable. 

One can show that the two eigenvalues of e~"''^ F{T) obey the relation A1A2 = 1, reflect- 
ing the incompressible piecewise Hamiltonian flow. Therefore the eigenvalues are complex 
conjugates in the stable regions and real in the unstable regions. At the boundaries both 
eigenvalues equal either +1 or —1, in which case the system oscillates with period T or 2r, 
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respectively. Therefore the trace of the Floquet operator is ±2 at the stabihty boundaries, 
leading to the following boundary equation: 

In the absence of damping the solutions depend only on A and the ratio r = uJo/tOp, 
and these particular results are shown by the dashed curves in the first panel of Fig. |l|. 
The regions of instability appear as "tongues" starting from integer and half-integer values 
of r at small modulation amplitudes. Note that the stability boundaries surrounding the 
tongues that emerge from integer values of r correspond to both Floquet eigenvalues being 
equal to 1, while those surrounding tongues that emerge from half-integer values of r are 
associated with the Floquet eigenvalues equal to —1. Note that no abrupt transition (nor a 
visible transition of any kind) is seen at the line A = 1 which marks the transition between 
an oscillating-oscillating behavior and a saddle-oscillating one. As one would expect, on the 
other hand, the system is nearly always unstable for sufficiently small ujp and sufficiently 
large A (upper right corner of the figure) . 

Acknowledgements: We would like to thank C. Van den Broeck and R. Kawai for stim- 
ulating discussions. This work was supported in part by the National Science Foundation 
under grant No. PHY-9970699. 

References 

[1] L. D. Landau and E. M. Lifschitz, Mechanics (Pergamon Press, Oxfors, 1969). 

[2] V. I. Arnold, Equations Dijferentielles Ordinaires (Editions Mir, Moscou, 1974); En- 
glish translation: Ordinary Differential Equations (Springer- Verlag, Berlin, 1992); 
Mathematical Theory of Classical Mechanics (Springer- Verlag, Berlin, 1989). 

[3] A. H. Nayfeh and D. T. Mook, Nonlinear Oscillations (Wiley, New York, 1979). 

[4] C. Van den Broeck and I. Bena, "Parametric resonance revisited," preprint. 

[5] E. Kh. Akhmedov, J. Physics 54, 47 (2000). 

[6] M. Calvo, Phys. Rev. B 60, 10953 (1999). 

[7] I. Zlatev, G. Huey and P. J. Steinhardt, Phys. Rev. D 57, 2152 (1998). 

[8] J. W. Miles, J. Fluid Mech. 148, 451 (1984); S. T. Milner, J. Fluid Mech. 225, 81 
(1991). 

[9] S. Riyopoulos, J. Plasma Phys. 46, 473 (1991); Phys. Rev. Lett. 68, 3303 (1992). 
[10] M. Calvo, Phys. Rev. B 55, 10571 (1997). 

[11] N. V. Fomin, O. L. Shalaev and D. V. Shantsev, J. Appl. Phys. 81, 8091 (1997); K. 
Otsuka, D. Pieroux, W. Ting-Yi and P. Mandel, Optics Lett. 22, 516 (1997). 



18 



[12] M. Schienbein and H. Gruler, Phys. Rev. E 56, 7116 (1997); F. S. Prato, M. Kavaliers, 
A. P. Cullen and A. W. Thomas, Bioelectromagnetics 18, 284 (1997). 

[13] D. Ballon et. al.. Magnetic Resonance in Medicine 39, 789 (1998). 



[14] R. Kobes and S. Peles, |nonl-phys 0005005| (2000). 

[15] I. Bena and C. Van den Broeck, Europhys. Lett. 48, 498 (1999). 

[16] J-Y. Ji and J. Hong, J. Phys. A 31, L689 (1998). 



19 



